Spin-Orbit Coupling and Ion Displacements in Multiferroic TbMnOs 
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The magnetic and ferroelectric (FE) properties of TbMnOs were investigated on the basis of 
relativistic density functional theory (DFT) calculations. We show that, due to spin-orbit coupling, 
the spin-spiral plane of TbMnOs can be either the be- or afe-plane, but not the ac-plane. As for 
the mechanism of FE polarization, our work reveals that the "pure electronic" model by Katsura, 
Nagaosa and Balatsky is inadequate in predicting the absolute direction of FE polarization. Our 
work indicates that to determine the magnitude and the absolute direction of FE polarization in 
spin-spiral states, it is crucial to consider the displacements of the ions from their centrosymmetric 
positions. 
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Recent studies on magnetic ferroelectric (FE) materi- 
als have shown that electric polarization can be signif- 
icantly modified by the application of a magnetic field 
P, i, i, i, i, i, 0]. Perovskite TbMnOa with a spin- 
spiral magnetic order is a prototypical multiferroic com- 
pound with a gigantic magnetoelectric effect Cur- 
rently, there are two important issues concerning the FE 
polarization of TbMnOa. One concerns the origin of FE 
polarization. Model Hamiltonians studies of spin-spiral 
multiferroic compounds have provided two different pic- 
tures. In the Katsura-Nagaosa-Balatsky (KNB) model 
, the hybridization of electronic states induced by spin- 
orbit coupling (SOC) leads to a FE polarization of the 
charge density distribution even if the ions are not dis- 
placed from their centrosymmetric positions. In contrast, 
the model study by Sergienko and Dagotto Q concluded 
that oxygen ion displacements from their centrosymmet- 
ric positions are essential for the FE polarization in mul- 
tiferroic compounds [1]. When carried out with the ions 
kept at their centrosymmetric positions, density func- 
tional theory (DFT) calculations [QJ for the spin-spiral 
states of LiCuV04 predict FE polarizations that agree 
reasonably well in magnitude with experiment [lo| , which 
is in apparent support of the KNB model It is, there- 
fore, important to check which model, the KNB or the 
"ion displacement" model, is relevant for the FE polariza- 
tion in TbMnOa. The other issue concerns the spin-spiral 
plane of TbMnOs. Under a magnetic field, the spin-spiral 
plane of TbMnOa can be either the &c-plane or the ab- 
plane, but not the ac-plane. To explain this observation, 
it is necessary to probe the magnetic anisotropy of the 
Mn'^+ ion. The magnetic anisotropy of the Tb^+ ion 
might be also relevant for the magnetoelectric effect, as 
suggested by Prokhnenko et al. 11]. 

In this Letter, we investigated these issues on the ba- 
sis of DFT calculations and found that the considera- 
tion of the ion displacements is essential for the FE po- 
larizations in the spin-spiral state of TbMnOa, and the 
KNB model can be erroneous even for predicting the ab- 



solute direction of FE polarization. The absence of the 
ac-plane spin-spiral in TbMnOs is explained by the mag- 
netic anisotropy of the Mn'^+ ion. 

Our calculations were based on DFT plus the on-site 
repulsion U metho d I l2l| within the generalized gradi- 
ent approximation il3|. We used U^ff = 6.0 eV on Tb 
4f states. With other Ueff values for Tb, similar re- 
sults were obtained. For Mn 3d states, we employed 
Ueff = 2.0 eV, which leads to the spin exchange in- 
teractions between the Mn'^+ ions that are consistent 
with the observed magnetic structure of TbMnOa (see 
below) . For the calculation of FE polarization, the Berry 
phase method [3] encoded in the Vienna ab initio sim- 
ulation package (VASP) was employed 15, l3|, in which 
the Tb 4f electrons were treated as core electrons. For the 
study of the SOC effect associated with the Tb 4f elec- 
trons, we used the full-potential augmented plane wave 
plus local orbital method as implemented in the WIEN2k 
code 17[. Due to the small value of the spin anisotropy, 
we employed the convergence threshold of 10~^ for elec- 
tron density. As shown in Fig. [2 the experimental crys- 
tal structure of TbMnOs has a distorted GdFeOa- 
type orthorhombic perovskite structure with space group 
Pbnm. In our calculations, the experimental structure 
was used unless otherwise stated. 

The basic electronic structure of TbMnOs is shown 
in Fig. [5] in terms of the density of states (DOS) of the 
FM state calculated by WIEN2k with the Tb 4f states 
treated as valence states. The system is an insulator with 
energy gap of about 0.5 eV between the spin-up Mn 3d 
eg states. The spin-up Tb 4f states are all occupied while 
the spin-down Tb 4f states occur at around —2.8 eV with 
the remaining down-spin Tb 4f states above the Fermi 
level. This feature is consistent with the configuration 
of the Tb'^"'' ion. Partial DOS analysis indicates that 
the occupied down-spin states have mainly the fz{x^-y^) 
orbital character with a slight contribution from the fxyz 
orbital. 

For the spin exchange between Mn'^^ ions, there are 
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FIG. 1: (Color online) Perspective view of the orthorhombic 
structure of TbMnOs. The large, medium, and small spheres 
represent the Tb, Mn, and O ions, respectively. The Mn-Mn 
spin exchange paths Jab, Jaa, Jbb, and Jcc are also indicated. 
The solid vectors denote the easy axes for the Tb^+ and Mn^+ 
spins. 
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FIG. 2: (Color online) Electronic structure obtained for the 
FM state of TbMnOs from the WIEN2k calculations. The 
total DOS, the Tb 4f partial DOS, and the Mn 3d partial 



DOS are shown, 
indicated. 
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four spin exchange paths to consider as shown in Fig. [T] 
Jab is the intralayer {ab plane) nearest neighbor (NN) ex- 
change interaction, Jaa and Jbb are the intralayer next- 
nearest neighbor (NNN) exchange interactions along the 
a and b directions, respectively, and Jcc is the interlayer 
NN exchange interaction along the c direction. The high- 
spin Mn^+ ions have the tig^g configuration with the Cg 
states forming a staggered orbital ordering of the c?3a;2_r2 
and d-^y2_j.2 states. Such an orbital ordering induces an 
NN intralayer ferromagnetic (FM) exchange, and an NN 
interlayer antiferromagnetic (AFM) exchange. The co- 
operative Jahn- Teller distortion associated with the or- 
bital ordering leads to a large NNN intralayer super- 
superexchange along the b direction [ll,[2^. Our VASP 
calculations lead to Jab = —1.52 meV, Jaa = 0.57 meV, 
Jbb — 0.85 meV, and Jcc = 0.50 meV. The calculated 
Jab, Jbb and Jcc values are consistent with those ex- 
pected for TbMnOa Jaa was predicted to be FM 



by Kimura et al. [l9|, but our calculations show it to 
be AFM. The classical spin analysis based on the Freiser 
method [2l| predicts that the spin ground state is a spin 
spiral with the modulation vector q = (0, gj,,0) with 
Qy — acos(— 5^)/7r = 0.15. The prediction of a spin- 
spiral ground state is in agreement with the experimental 
observation p^ . which shows an incommensurate spiral 
with q = (0, 0.27, 0) below 28 K. It is worthwhile to point 
out that if a smaller Uef / is used Jbb will be stronger, and 
Jab be weaker, hence leading to a larger qy. 

To examine the magnetic anisotropy of the high-spin 
Mn'^+ ion in TbMnOs, it is necessary to consider DFT-I-U 
calculations with SOC included. TbMnOs has four Mn^+ 
ions per unit cell. The easy axes of these ions are not 
coUinear but related by symmetry. Since the effect of 
SOC is largely local in nature, we consider for simplic- 
ity only one Mn ion (i.e., Mnl in Fig. [T]) per unit cell 
by replacing the remaining three Mn^+ ions with Sc'^^ 
ions that have no magnetic moment. To remove any pos- 
sible coupling with the Tb 4f moment, we replace the 
Tb'^+ ions of the unit cell with La^^^ ions that have no 
magnetic moment. The energy dependence upon the Mn 
spin direction obtained from DFT-f U-I-SOC calculations 
is shown in Fig. ^a) for the cases when the Mn spin lies 
in the ab and be planes. The energy minimum occurs at 
60° and 70° for the case when the spin lies in the ab and 
be planes, respectively, and the energy minimum for the 
a6-plane is lower than that for 6c-plane by 0.12 meV/Mn. 

The above results can be readily understood by ana- 
lyzing the effect of the SOC Hamiltonian [Ulijl. As 
shown in Fig. [SJb), each MnOe octahedron of TbMnOa 
is axially elongated. With one of the two longest Mn-0 
bonds taken as the local z axis and neglecting the slight 
difference between the other four short Mn-0 bonds, the 
d-block levels of the Mn'^^ (d^) ion are described as 
shown on the right hand side of Fig.[3ljb), which leads to 
the electron configuration {xzY {yzY {xyy {z'^Y {x"^ —y"^)^ . 
The SOC Hamiltonian AL • S leads to an interaction of 
the empty d^2_y2 (eg) state with the other d-states. The 
strongest SOC occurs between two Cg states {dxy and 
dr^2_y2) with the maximum eiiergy lowering when the 
spin lies in the local z direction [2J] . For the Mnl ion, the 
6 and </) angles of the local z direction in the global coor- 
dinate system are 9 = 80.2° and (/> = 60.5°, respectively. 
In good agreement with these values, our SOC analysis 
based on tight-binding calculations 25| shows that the 
actual easy axis of Mnl has = 84° and (p = 60°. Thus 
the easy axis is close to the ab plane, and is far from the 
c axis. This explains why the energy minimum for the 
Mn spin lying in the a6-plane has a lower energy than 
that lying in the 6c-plane. In addition, the easy axis is 
closer to the b axis than to the a axis. Our result is con- 
sistent with the experimental observation [2^ that the 
sine- wave modulation of the Mn magnetic moment below 
40 K has the direction parallel to the 6-axis. Since the 
easy axis of the Mn spin is far from the c direction, it is 
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FIG. 3: (Color online) (a) The dependence of energy on the 
direction of the Mn spin in the afe-plane or the bc-plane. (b) 
The local structure of the distorted MnOe octahedron. The 
numbers give the Mn-O bond lengths in A. The local coor- 
dinate system is also indicated. The right-hand-side diagram 
illustrates the electron configuration of the Mn'^^ (d**) ion. 
The label "SOC" denotes that the largest SOC mixing occurs 
between the d^y and d^i_yi states, (c) The dependence of 
energy on the direction of Tb spins in the afe-plane. 



expected that the spin-spiral plane of TbMnOa is either 
the a5-plane or the 6c-plane, but not the ac-plane 
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Experimentally [24], it was found that the Mn moments 
of TbMnOs form a 6c-plane elliptical spiral with mb(Mn) 
= 3.9 > mc(Mn) = 2.8 fiB- This can be readily ex- 
plained by the fact that the easy axis of the Mn spins is 
close to the b axis, but far from the c direction. 

The magnetic anisotropy of the Tb'^^ ion in TbMnOa 
was also calculated in a similar manner; three of the four 
Tb^"*" ions in a unit cell were replaced by La''+ ions with 
all the Mn'^^ ions of a unit cell replaced by Sc'^^ ions. 
Our DFT-fU+SOC calculations show that the Tb spin 
prefers to lie in the ab plane since the state for the spin 
parallel to the c-axis is higher than that for the spin in the 
ab plane by at least 0.83 meV/Tb. The energy calculated 
for Tbl (as labeled in Fig. [T] as a function of the angle 
(j) with = 90° is presented in Fig. [3]Jc), which reveals 
that the energy miminum occurs at around 4> = 145°. 
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FIG. 4: (Color online) The atomic displacements (indicated 
by solid arrows) of the spin-spiral states after geometry opti- 
mization for (a) the fee-plane spiral and (b) the afe-plane spiral. 
The dashed arrows indicate the directions of the Mn spins in 
the spiral states. The directions of the electric polarizations 
are also shown. 



The large anisotropy energy is consistent with the Ising 
behavior of Tb moment . The easy axes of the Tb"^"*" 
ions presented in Fig. [T] show that they lie symmetrically 
around the 6- axis with the angle 55°. This is in excellent 
agreement with the observed angle of 57° [2^ . 

The electric polarization of TbMnOa in the spin-spiral 
state was calculated using the VASP. To reduce the com- 
putational task, we considered the q = (0, 1/3,0) state, 
which we simulated by using a 1 x 3 x 1 supercell. Using 
the experimental centrosymmetric structure, the electric 
polarization from the pure electronic effect was calcu- 
lated to be P = (0,0,0.5) \iClm^ [H, S [s^ for the 
&c-plane spiral shown in Fig. (Ha). Experimentally, the 
magnitude of the electric polarization for the 6c-plane 
spiral is about 600 /iC/771^ [3l[, which is three orders 
of magnitude larger than the value calculated with no 
geometry relaxation. For the a6-plane spiral shown in 
Fig. mb), the electric polarization is calculated to be 
P = (—331.0,0,0) [iCIrr? in the absence of geometry 
relaxation. The absolute directions of the FE polariza- 
tions obtained for the be- and a6-plane spin-spiral states 
from DFT calculations without geometry relaxation are 
opposite to those predicted by the KNB model (Eq. 6 of 
ref. Q). To examine the effect of ion displacements in the 
spin-spiral states on the FE polarization, we optimized 
the atom positions of TbMnOa in the be- and a5-plane 
spiral states by performing DFT+U-f SOC calculations 
and calculated the electric polarizations of TbMnOa us- 
ing the relaxed structures. These calculations lead to Pc 
= —424.0 [iC jrr? for the 6c-plane spiral, and to Pa = 
131.2 fj,C/m^ for the a6-plane spiral. The absolute di- 
rections of the FE polarizations are switched for both 
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spin-spiral states under geometry relaxation. The cal- 
culated FE polarization for the 6c-plane spiral using the 
relaxed structure is now much closer to the experimental 
value. Furthermore, the direction of the FE polarization 
is in agreement with experiment 28, 2^ 3^, 32 1 as well 
as the KNB model. However, this agreement between 
the KNB model and experimental values is fortuitous; 
for LiCuV04 and LiCu2 02, which consists of Cu02 rib- 
bon chains, the KNB model predicts the wrong absolute 
direction of the FE polarization for the a6-plane spiral, 
but the correct absolute direction of the FE polarization 
for the 6c-plane spiral, regardless of whether the unre- 
laxed or the relaxed crystal structures are employed 
The failure of the KNB model could be due to the fact 
that it was derived for the t^g systems while the tg states 
are important in Mn'^+ and Cu^+ systems. Furthermore, 
we notice that the FE polarization of the a6-plane spi- 
ral calculated by using the relaxed structure, now along 
the a direction, is smaller in magnitude than that calcu- 
lated by using the experimental structure. In the case of 
LiCuV04 and LiCu202, the geometry relaxation in the 
spin-spiral states was found to enhance the magnitudes 
of the FE polarizations without changing their absolute 
directions [9|. 

To show how the ion displacements break the cen- 
trosymmetry, we present the pattern of the atom dis- 
placements with respect to the centrosymmetric struc- 
ture in Fig. m For the 5c-plane spiral, the O atoms at 
the Wyckoff 4c position (i.e., 01) lying in the ab plane 
of the Tb atoms almost do not move. The movement 
of all the O atoms at the Wyckoff 8d position (i.e., 02) 
close to the ah plane of the Mn atoms have a component 
along the c direction. All the Mn and Tb atoms have 
a displacement along the — c direction. The largest dis- 
placements along the c direction (about 2.4 x 10~^ A) 
occur at the Mn sites rather than at the O sites. This 
finding is in contradiction to the assumption introduced 
in the "ion-displacement" model Considering that 
the Born effective charge is positive for the Tb'^+ and 
Mn'^+ cations, and negative for the 0^~ anion, the total 
electric polarization is expected to be along the — c di- 
rection, in agreement with the DFT calculation. For the 
a6-plane spiral, some O atoms have displacements along 
the —a direction, but other O atoms have displacements 
along the a direction. The sum of the O ion displace- 
ments is along the —a direction. The occurrence of alter- 
nating O displacements is consistent with the prediction 
made by Sergienko and Dagotto [3], and is responsible 
for the smaller electric polarization when compared with 
the case of the 6c-plane spiral. All the Mn atoms have a 
displacement along the a direction. Another unexpected 
finding is that some Tb atoms have the largest displace- 
ments with a large component along the h or —6 direction 
and a small component along the a direction. 

In summary, the absence of the ac-plane spin-spiral in 
TbMnOa is explained by the magnetic anisotropy of the 



Mn^+ ion. The calculated easy axis for Tb is in excellent 
agreement with the experimental result. The considera- 
tion of the ion displacements in the spin-spiral states of 
TbMnOa is essential in determining the magnitude and 
the absolute direction of the FE polarizations, which is in 
support of the "ion-displacement" model. Surprisingly, 
the displacements of the Mn'^+ and Tb'^+ ions are gener- 
ally greater than those of the 0^~ ions. The KNB model, 
however, can fail to describe both the magnitude and the 
absolute direction of FE polarization. 
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